The First Law for Slowly Evolving Horizons 
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We study the mechanics of Hayward's trapping horizons, taking isolated horizons as equilibrium 
states. Zeroth and second laws of dynamic horizon mechanics come from the isolated and trapping 
horizon formalisms respectively. We derive a dynamical first law by introducing a new perturbative 
formulation for dynamic horizons in which "slowly evolving" trapping horizons may be viewed as 
perturbatively non-isolated. 



The laws of black hole mechanics are one of the most 
remarkable results to emerge from classical general rel- 
ativity. Recently they have been generalized to locally 
defined isolated horizons 0. However, since no matter 
or radiation can cross an isolated horizon, the first and 
second laws cannot be treated in full generality. Instead, 
the first law arises as a relation on the phase space of 
isolated horizons rather than as a truly dynamical rela- 
tionship. Even existing physical process versions of the 
first law [HQ consider transitions between infinitesimally 
separated isolated (or Killing) horizons. In this paper, 
we introduce a framework which allows us to extend the 
first law to all slowly evolving horizons, even those whose 
initial and final states are not infinitesimally separated. 
Additionally, we provide a simple characterization of how 
close a horizon is to equilibrium, which will be useful in 
the final stages of numerical simulations of black hole 
collisions. To obtain such a law, we first need a local dy- 
namical definition of a black hole horizon. Hayward has 
introduced future outer trapping horizons (FOTHs) for 
exactly this purpose. Furthermore, he has shown that 
they necessarily satisfy the second law — their area can- 
not decrease in time 

In this paper, we will show that both the zeroth and 
first laws are also applicable to FOTHs. In order to do 
so, we introduce dynamical notions of surface gravity and 
angular momentum which are applicable to all such hori- 
zons. It follows immediately that the surface gravity is 
necessarily constant if the horizon is in equilibrium (i.e. 
isolated) . Next, we introduce the notion of a slowly evolv- 
ing horizon, for which the gravitational and matter fields 
are slowly changing. It is only in this limited context 
that we expect to obtain the dynamical first law. We 
will show that this is indeed the case. 

Let us begin by recalling Hayward's definition: 

Future Outer Trapping Horizon: A future outer 
trapping horizon (FOTH) is a smooth three-dimensional 
sub-manifold H of space-time which is foliated by a pre- 
ferred family of space-like two-sphere cross-sections H v , 
with future-directed null normals i and n. The expan- 



sion 9(g) of the null normal I vanishes. Further, both the 
expansion 9i n \ of n and C n Q(t) are negative. 

This definition captures the local conditions by which 
we would hope to distinguish a black hole. Specifically, 
the expansion of the outgoing light rays is zero on the 
horizon, positive outside and negative inside. Addition- 
ally, the ingoing null rays are converging. As we will see, 
these horizons can be space-like or null and include both 
equilibrium and non-equilibrium states — an important 
feature for our perturbative study. By comparison, if one 
is interested only in the dynamical phase, the space-like 
dynamical horizons recently introduced by Ashtekar and 
Krishnan |f| are more relevant. However, such horizons 
are always expanding and so not so suitable for studying 
transitions to and from equilibrium. 

Hayward has extensively studied the properties of 
FOTHs. Here, we summarize only those of his results 
which arc important for our work. First, consider the 
quantities associated with the null vector fields £ and 
n, which we normalize so that £ ■ n = — 1. We denote 
their relative expansions 9^ and 6i n y Their twists are 
zero since they are normal to the H v cross-sections of the 
horizon. Finally, we write their shears as and cr^ ■ 

Next, for each choice of the fields I and n, there exists 
a scalar field C on H so that 



V° = t - Cn a and r a = l a + Cn a , 



(1) 



are respectively tangent and normal to the horizon. Note 
that V • V = — r • r = 2C. Hayward |^| has shown that if 
the null energy condition holds, then 



C > 



(2) 



on a FOTH. Thus, the horizon must be either space-like 
or null, and the second law of trapping horizon mechanics 
follows quite easily. If q a b is the two-metric on the cross- 
sections, and y/q is the corresponding area element, then 



£vV9 = -C9( r . 



(3) 
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By definition 6*(„) is negative and we have just seen that 
C is non-negative. Hence we obtain the local second law: 
If the null energy condition holds, then the area element 
y/q of a FOTH is non-decreasing along the horizon 
Clearly on integration, the same law applies to the total 
area of the horizon. It is non-decreasing and remains 
constant if and only if the horizon is null. 

To further restrict the rescaling freedom of the null 
vectors, we require that C yv = 1, where v is the foliation 
parameter labeling the cross-sections. Thus, choosing the 
foliation parameter fixes the length of V a , r°, and the 
null vectors. However, at this stage we are still free to 
change the foliation labeling and, by doing so, rescale V a 
and the null vectors by functions of v. From the isolated 
horizon perspective, the normalization we have chosen is 
very natural, as in that case V a — £ a and it is customary 
to choose the foliation parameter v, so that n — —dv. 
Physical characterization of trapping horizons — The 
laws of black hole mechanics are given in terms of phys- 
ical quantities such as energy, angular momentum, sur- 
face gravity and area. One of the advances of the isolated 
horizon formalism was to provide definitions of all these 
quantities at the horizon, without reference to space-like 
infinity or the space-time in which the horizon is embed- 
ded In generalizing these definitions to FOTHs, we 
will be motivated by two requirements: i) the new defini- 
tions should match the old ones when trapping horizons 
are isolated and ii) the new definitions should depend 
only on quantities that are intrinsic to the horizon, and 
as such truly be properties of the horizon itself. That 
said, the ultimate justification for these expressions will 
be found in their utility in the calculations that follow. 

For isolated horizons both surface gravity and angular 
momentum were defined with respect to a one-form uj a 
which can be written as: 

u a := -n b V_ a £ b , (4) 

where the arrow signifies that the derivative is pulled 
back to the horizon. Then, the isolated horizon surface 
gravity is given by Kg = £ a uj a while angular momentum 
information is contained in the other components of u> a . 

Now, uj a written in this form (equation continues 
to be an intrinsically well-defined quantity for trapping 
horizons. Then, an obvious generalization is to define: 

k v := V a Lu a (= -n b V a V a V h ) and 

J v ■= — I d 2 xy^(p a Q a , (5) 

where tp a is any vector field tangent to the horizon cross- 
sections, and uj a is the projection of uj a into the 2-surface 
H v . In the case where ip a is divergence free, J v is in- 
dependent of the normalization of £ and equal to other 
popular expressions for angular momentum (such as the 
Komar, Brown- York 6j or Ashtekar-Krishnan || defini- 
tions). It is clear that both the surface gravity k v and 



angular momentum J v expressions reduce to their iso- 
lated horizon values if the horizon is null. 
The constraint law — Consider the set of vector fields 
in H that generate onc-paramctcr families of diffcomor- 
phisms that map two-dimensional cross-sections H v of 
the horizon into each other. Any such vector field may 
be written in the form X a = x V a + x a , for some func- 
tion x (v) and vector field x a that is everywhere tangent 
to an appropriate cross-section. Then, by integrating the 
r a X b component of the Einstein equations over H v , we 
obtain the following relationship: 




A full derivation of this result will be given in [3 and in 
another paper we will see that this relation plays a crucial 
role in the Hamiltonian formulation of general relativity 
on manifolds with FOTHs as boundaries ||. Here, we 
show that in certain restricted situations equation © 
becomes a first law for dynamical black holes. 
Quasi-stationary horizons — At first glance, one might 
think that © is already the dynamical first law of black 
hole mechanics. After all, it relates rates of change of 
area and angular momentum to fluxes across the horizon. 
However, there are several reasons why this is not so. The 
first and most important is that we should not expect the 
standard kcl form of the first law to hold in all dynamical 
situations. In thermodynamics, it is only in the quasi- 
static case that it is possible to write the energy balance 
equation as 

dE = TdS + work terms . (7) 

Furthermore, in the general case, there is no clear inter- 
pretation of the right hand side of © as a flux of energy 
through the horizon. Instead, we will require the horizon 
to be "quasi-stationary" and then obtain a first law. 

Heuristically, it is clear that the properties of a quasi- 
stationary horizon, such as the area and surface gravity , 
should be slowly varying. However, since there is a rescal- 
ing freedom in the vector field V, requiring (C\>\/q)/^/q 
to be small is not a meaningful condition — it can be 
satisfied on any trapping horizon by suitably rescaling V. 
Furthermore, we cannot fix the norm of V (or its average) 
to unity since we are interested in the limit as V becomes 
null. Instead, we introduce the one-form \ a := d a v satis- 
fying V • x = 1 , and require (£ y \fq) (C x Vg) /q to be small 
compared to the characteristic scale of the horizon. This 
condition is invariant under rescalings and the expression 
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vanishes whenever the horizon is null. On a space-like 
horizon, (£ v^q)(£ x *Jq) = qCO?-. and so is closely re- 
lated to the expansion of the surface (see equation J3J). 
By evaluating this expression, we can invariantly deter- 
mine the expansion rate of the horizon and choose the 
foliation accordingly. Then, we require other fields to be 
slowly varying with respect to this foliation. 
Slowly Evolving Horizon: A region of a future outer 
trapping horizon H with v G [fi,i>2], is slowly evolving 
(at a rate e) if there exists a parameter e <C 1 such that, 

1. On every cross section H v with v € [vi,?^], 
/ d 2 xj~q (q ab V a Vt) (q cd V c X d) < e 2 . 

2. The foliation parameter v is chosen so that |V| = 
\J2C ~ e, and its rates of change are similarly small 
over the horizon. 

3. The one-form Q a and expansion n are slowly evolv- 
ing: |£ v ^a| < e/r% and \C V 9^\ < e/r%, where 
Th is the area radius of the horizon, au = ^r 2 H . 



4. \n\ 



I 2 , |cr (n) | 2 and T ab n a n b 



l/rfi or smaller. 



The first condition gives an invariant characterization 
of slow expansion. The second condition restricts the fo- 
liation to guarantee that the area is slowly evolving with 
respect to V a , specifically (£yr#) ~ e 2 - The third re- 
quires other geometric quantities to also be slowly evolv- 
ing, while the fourth fixes a reasonable horizon geometry 
and demands that conditions in the surrounding space- 
time not be too extreme, here 1Z is the Ricci scalar of 
the two-surfaces H v . An immediate consequence of our 
definition is that isolated horizons are examples of slowly 
evolving horizons with e = C = (provided condition 4 
is satisfied). 

Let us now consider the implications of this definition. 
Here, we will simply state the consequences, a more com- 
plete discussion will be given in [7j. For concreteness, we 
restrict the allowed matter to be scalar and/or electro- 
magnetic fields. From a projection of the Einstein equa- 
tions and condition 2, we can show that 



M\ 



e and T ab £ a £ b 



(8) 



(Here and in the following, we choose to focus on the e 
factors. The powers of rn required to make the equa- 
tions dimcnsionally correct are omitted.) Then, the al- 



lowed form of the matter fields forces T ab 



and 



£v(T a b£ a n ) ~ e. Further application of the Einstein 
equations and Bianchi identities gives the Weyl compo- 
nents ~ e and CyTZ ~ e. On an isolated horizon 
each of these quantities would be zero. Lastly, condition 
3 is sufficient to guarantee that 



|£v««|~e and \q a "V& K v \ ~ e . 



(9) 



Therefore, on a slowly evolving horizon, the surface grav- 
ity is a slowly varying function and we can write 



..(()) 



6K 



(1) 



(10) 



where «(°) is a constant. Recall that an isolated horizon 
corresponds to e = 0. Hence, we immediately obtain the 
zeroth law: the surface gravity of a FOTH is constant if 
the horizon is isolated. 

The dynamical first law — We can now derive the first 
law for slowly evolving horizons. First we consider evolu- 
tion with respect to V a , an appropriate evolution vector 
for a horizon with no obvious axis of symmetry — a non- 
rotating horizon. Setting x a = and x Q — 1 in the 
constraint law 10 and expanding in powers of e, we find 
that all terms vanish at zeroth and first order. Then, to 
0(e 2 ), the first law for a slowly evolving horizon reads: 



1 



8ttG 



E 



/ d 2 x^f\ 



T ab e a e b + 



1 



8ttG' 



(11) 

where the dot signifies a derivative with respect to n. It 
says that the surface gravity multiplied by the change in 
area is equal to a flux of energy through the horizon. This 
flux is comprised of two terms, both of which are positive. 
The first is the flux of matter through the horizon and 
is familiar from standard physical process versions of the 
first law The second term is a flux of gravitational 
shear through the horizon — which would naturally be 
interpreted as a flux of gravitational radiation. A similar 
term has been obtained previously by Hawking 9] when 
considering perturbations of an event horizon. 

We would like to "integrate" (| 1 1 11 in order to obtain an 
expression for the energy of the horizon. To this end, we 
further restrict the choice of V a by requiring that: 



1 

27^" 



(12) 



Note that while we cannot prove that this rescaling is 
possible on every slowly evolving horizon, we can show 
that if the horizon satisfies a certain genericity condition, 
it will be. In particular, perturbations to Schwarzschild 
and non-extremal Kerr horizons will satisfy the condi- 
tion. Since is a function of ru alone, we can inte- 
grate the k a term in i|ll|) to obtain the usual expression 
for the energy 



2G 



(13) 



Rotating horizons — Next, consider rotating horizons. 
In (0) we gave a definition for the angular momentum. 
However, this quantity only has a physical significance if 
the vector field ip a is an (approximate) symmetry of the 
horizon. Hence, we define an 

Approximately symmetric horizon: The vector field 
ip a is an approximate symmetry of a trapping horizon if 
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1. (f a is tangent to the cross sections of the horizon, 
generates a family of closed integral curves and is 
normalized so that those curves have afhne param- 
eter length 2ir. 

2. tp a is Lie dragged up the horizon by V a , i.e. 
£v<P a = 0. 



3. tp a is an approximate symmetry of the horizon 
geometry and matter fields: \q ab V a fb\ < e 2 , 

< e and 



v q a b\ < e, 



< e, 



<P °(n) I 



< e 



(These conditions are analogous to 



those satisfied by V a on a slowly evolving horizon) . 

For a slowly evolving horizon which is approximately 
symmetric, the angular momentum J v is a meaningful 
quantity. Going back to the constraint law and setting 
x a = and x a — ip a we find its time rate of change. As in 
the non-rotating case, equation © will vanish at zeroth 
and first orders in e. To second order we have: 



H, 



T^o. , b 
abT <fi 



1 



(14) 

Thus, the rate of change of angular momentum of the 
horizon depends upon the flux of matter and gravita- 
tional fields through the horizon. 

To get a general first law, we would like to combine 
the rate of change of angular momentum with the first 
law for non-rotating horizons l|f f|) . To do this, we again 
fix the average value of the surface gravity as well as an 
angular velocity Q (we are following a similar strategy to 
that taken in though with the caveat that here only 
the average value of the surface gravity is fixed). We do 
this by requiring that both take the same values as in a 
Kerr space-time with the same and an'- 



ft := 



2GJ E 



and 



' H 



4G 2 J 2 



H 



2r| 



H 



(15) 

With area and angular momentum slowly changing up 
the horizon, the angular velocity is also slowly varying. 

Finally, we consider the full constraint law © with the 
evolution vector field X a chosen as 



of the equation to obtain as an expression for the energy 



E = 



H 



AG 2 J 2 r 



2Gr 



H 



(18) 



This E is equal to the energy of a rotating isolated hori- 
zon with the preferred choice of normalization for £, and 
in particular is the mass of a Kerr black hole with pa- 
rameters an and J v . 

The notion of a slowly evolving horizon provides a 
bridge between the equilibrium of isolated horizons and 
the fully dynamical horizons of . In particular, our first 
law HI7J) bears a striking resemblance to the dynamical 
horizon energy balance formula. There are, however, sev- 
eral important differences. Firstly, we restrict to slowly 
evolving horizons, which are near to equilibrium. Doing 
so enables us to provide a locally defined surface gravity 
k v which is shown to be slowly varying over the horizon. 
This surface gravity maintains its usual interpretation 
as the acceleration of a vector field, V a , along the hori- 
zon. Secondly, we do not obtain an equivalent of the |C| 2 
term of [f| , which is likely related to angular momentum. 
However, we expect that term would vanish at the order 
of perturbation theory which we are considering. In the 
future, we plan to examine the connection between these 
two formulations in more detail and determine in what 
precise sense our first law can be seen as a perturbative 
form of the dynamical horizon energy balance formula. 

To summarize, we have examined conditions for which 
a slowly evolving FOTH may be said to be "perturba- 
tively non-isolated" . This has allowed us to obtain a truly 
dynamical version of the first law of black hole mechanics. 
However, we expect that the notion of a slowly evolving 
horizon will also find application in numerical studies of 
how horizons settle down to equilibrium. Specifically, one 
could track the approach of the parameter e to zero. 
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t a := V a + IV 1 ■ 



Then, we can once again expand out the constraint law 
© to order e 2 and this time obtain: 



E 
E 



InG 



K^°'d + nj v where 



d x\l q 



T ab t a r b 



1 



I6ttG 



r (JL-tQab) 



Since and Q are specific functions of only rjj and J v , 
it is once again possible to "integrate" the left hand side 
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